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Summary 

An experimental method for measuring rolling friction 
is proposed. The method is mechanically simple. It is 
based on an oscillator in a uniform magnetic field and 
does not involve any mechanical forces except for the 
measured friction. The measured pickup voltage is 
Fourier analyzed and yields the friction spectral response. 
The proposed experiment is not tailored for a particular 
case. Instead, various modes of operation, suitable to 
different experimental conditions, are discussed. 


Introduction 

Rolling friction measurements usually involve 
elaborate mechanical arrangements and, in some cases, 
are not easy to make (refs. 1 to 4). Moreover, there is a 
suspicion that dissipative or nondissipative forces, apart 
from the friction to be measured, affect the results. In 
this report we are proposing a new method of rolling 
friction measurement that does not rely on mechanical 
devices. The only moving part is the rolling body whose 
friction (against a plane surface) is to be measured. It 
involves the use of an oscillator and the measurement of 
magnetic forces. Both the analysis of the concept and the 
experimental approach are given. 


Analysis 

Mechanically the apparatus is extremely simple, it 
includes a rolling oscillator over a horizontal plane. The 
restoring torque t is obtained by a magnetic field B 
interacting with a magnetic dipole /T according to the 
following equation (ref. 5): 

T=ii'xB ( 1 ) 

The magnetic dipole results from dc current circulating in 
a coil placed inside the oscillating body. The directions of 
^and ^are parallel at equilibrium. If we assume small- 
angle oscillations, the equation of motion for the angle 
0(t) is 

le + BiiO+f^O (2) 

where I is the moment of inertia and / is the rolling 
friction torque. 

In cases of an electrically conducting oscillator, losses 
due to induced currents may be important. These losses, 
which are not included in equation (2), can be sub- 
stantially reduced by reducing the conducting portions of 
the oscillating body. If the material whose friction is to be 
measured is metallic, it can be attached as a thin metallic 


foil on the surface of a bulk-insulating oscillator. 
Bringing the induced currents to aminimum is important 
for minimizing the experimental error. If that cannot be 
done directly, it can be done electrically by means of a 
feedback circuit, as discussed later under experimentEil 
considerations. Although it is possible to practically 
eliminate induced currents, we shall include their effect in 
our preliminary analysis. Their net effect is a magnetic 
moment normal to ^and proportional to B ^/p, where p 
is the electrical resistivity. Thus the induced current term 
is a “Newtonian viscous damping” term. With it the 
equation of motion becomes 

ie+-B^+fiBe+f=o (3) 

P 

where a is a constant dependent on the geometry. With 
/=0, the “frictionless” or “pure viscous damping” case, 
the solution of equation (3) expresses exponentially 
decaying oscillations. When the friction is included, the 
equation is not soluble in the general case. Den Hartog 
(ref. 6) first solved accurately the externally driven 
(forced) oscillator with Coulomb (square wave) friction 
and combined Coulomb and viscous friction. 

Equation (2) or (3) expresses / as a differential 
equation of 6. The angle d{t) could be measured in 
different ways. We suggest the use of a magnetic pickup 
as the most suitable for our experimental approach. As is 
depicted in figure 1 a pickup coil Lp is placed in the 
oscillator so that its axis is normal to ^at equilibrium. 
The coils L, Lp, and are the restoring, pickup, and 
deriving coils, respectively. 

The voltage V(t) of the pickup coil is 

V{t) = - ^NpApB cos e e (4) 

where Np and Ap are, respectively, the number of turns 
and the cross-sectional area of Lp, For small-angle 
oscillations we can approximate 

V(t) = - \NpApBe (5) 

and 

ao\* V{T)dT¥aiV(t)+a2V{t)+f=0 (6) 

Jo 

where aj =0 with no induced currents. 

We would like to stress that although the proposed 
experiment is introduced herein with specific features, the 
frictional oscillator is far more general. Thus, for 
example, the restoring force could be gravity or another 
force, and d(t) could be measured by optical or 




mechanical means. Most of the discussion and analysis 
presented herein could be applied to any dissipative 
oscillator. Even the rolling friction is not essential and the 
method, in principle, is applicable to sliding friction as 
well. A specific apparatus is proposed because we believe 
it is the simplest way to apply the discussed method. 

Even though equation (3) (or eq. (6)) is soluble in 
certain cases, it is not soluble in general. Rather than 
approximating it by using a model for f, we invert the 
equation of motion (i.e., use the measured velocity d(t) 
to numerically obtain the friction /). This function /(/) 
can in principle be interpreted in terms of features 
inherent to the particular surfaces under study, ignoring 
features resulting from the main resonance If the 

coefficients a/ are known together with V(t),f(t) can be 
found directly (by using numerical differentiation and 
integration). We believe that this straightforward 
approach, even though right in principle, does not yield 
the greatest accuracy. This is so because numerical 
differentiation usually yields rather “noisy” results 
unless incorporated with smoothening routines. The 
smoothening, however, is limited to the extent that it 
does not smear out real features. The alternative method 
of Fourier (or Laplace) transforms used herein is 
potentially better in terms of discrimination between 
signal and noise. Also, it provides a way of eliminating 
aj, which is otherwise the main cause for systematic 
error. This is so because, as f, it represents a loss and 
affects the oscillatory decay. 

The coefficient a\ in equation (6) should be referred 
to as an experimentally adjustable parameter. The 
adjustment can be done by means of feedback electronic 
circuitry, which is described in the next section. If ai is 
adjustable, we would like to consider in detail two cases 
of special interest. One case is when ai=0, which 
describes either the zero induced currents or the induced 
currents compensated by the feedback system. The other 
case of special interest is when the energy fed back 
compensates for both the induced current loss and the 
friction loss. 


The Laplace transformation of equation (6) yields 
(ref. 7) 

ctViO =aVo(t) [l - (7) 

where a^(/) = [ y(t) dt and V(j{t) is the friction- 

less solution for V{t) (what V{t) would have been if 
/=0). In equation (7), t = 0 corresponds to the meiximum 
pickup voltage The frictionless Ko(0 for nj =0 is an 
undamped harmonic oscillator. Its determination 
involves two parameters, the oscillation period and 
amplitude, which are well approximated by the measured 
period and first-cycle amplitude of V{t) . The constant Q 2 
should be determined independently. Thus for zero 
induced currents the experimental procedure involves 
exciting the oscillator by a pulse in measuring V(t), 
extracting from it and olVU), and finally obtain- 

ing/by equation (7). 

If viscous losses exist {a^ 5*^0), determining V^it) may 
be a problem, since the oscillation decay rate is affected 
by both / and a\. It is always advantageous to minimize 
induced currents by reducing the metallic parts of the 
oscillator and enhancing the resistance of the current 
loops. In addition, we can compensate for these losses by 
feedback to We notice that the poles of aV^it) must 
be poles of aV{t) . Thus the feedback should be adjusted 
to result in an imaginary pole of aVQ (and aV). 
Practically this implies adjusting the feedback to get the 
best obtainable maximum of the half-line Fourier 

transform ( K(r) dr. To avoid being in the 
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vicinity of the aV poles whose inherent source is af 
(rather than aFo), check for the maximum close to the 
expected Kq resonance, namely cjq = We 

conclude that it is feasible to work in the zero a\ mode in 
any practical case and to obtain the half-line Fourier 
transform of /, 

It is quite interesting to deduce from equation (7) and a 
known theorem of Laplace transforms that (ref. 7) 

V{t) = Voit) - (V(r) Voit-r) dr (8) 

Thus the deviation is proportional to the 

convolution of the friction / with the frictionless function 
Fq. This conclusion, even though not useful for the 
present analysis, improves our understanding. It can be 
generalized to nonoscillatory problems. It is valid for any 
frictional problem whose frictionless part can be written 
as a linear and homogeneous differential equation. 

So far we have dealt with the decaying oscillation mode 
of operation. It is worthwhile to consider modes of 
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operation that maintain steady-state oscillations by 
providing the energy losses from an external source. At 
first we assume that the energy return is provided by a 
feedback force proportional to F, namely by adjusting a\ 
in equation (6) to maintain steady-state oscillations. The 
periodicity suggests the use of the finite Fourier 

transform £■„( F) = \ e'«"^F(T)</(cor). Transforming 

-X 

equation (6) we get 

aoE„ ( V) + inuHZiEn ( V) 

— a 2 ri^(t)^E„ (V) + inu}E„ (/) = 0 (9) 

By writing separate equations for the real and imaginary 
parts and denoting A„ = (ReE„(V), B„ = 9mE„(V), 
C„ = <ReE„(f), D„ = 9mE„(f), we have 

CO? — n2o)2 

C„= -aiA„-a 2 -^ B„ (10a) 

ncx) 


resonances (or different features) inherent to the 
frictional response. Typical relaxation times associated 
with the response of the particular surfaces in study 
should be reflected in £„(/). Further insight can be 
achieved by varying wq (by means of the restoring force 
dc current) and obtaining £*„(/, wq). 

Finally, we would like to point out the possibility of 
externed excitation with fixed amplitude and frequency 
(i.e., the “forced oscillator”). The equation of motion 
(without induced currents) is 

ao [ V t/r + 02 F+/= Q cos Qt (12) 

Jo 

which transforms to 


Dn= -aiBn + 02 


CO? — 


"0 


nco 


(10b) 


where n= 1, 3, 5, . . . Even /z’s are excluded because of 
the problem symmetry, since V(t) is antisymmetric with 
respect to a translation by a half-oscillation period. 
Equations (10) provide the spectral distribution of f 
versus the measured function F(r) and the parameters 
wq, 02^ Sind Q\. For zero induced currents a\ can be 
obtained by measuring the feedback current to For 
substantial induced currents, however, an accurate 
determination of a\ is hard and this mode of operation is 
undesirable. 

If induced currents do not exist, we can simplify the 
solution by filtering the feedback current and keeping 
only its first harmonic. Then a\V in equation (6) is 
replaced by b sin o>t, and we get for odd n 


Cn=-02— 


no) 


(11a) 


CO? — /l2co2 

— mo ^ ^ 


(11b) 


where j is the Kronecker function. Similarly to the case 
of decaying oscillations, equation (11) suggests that 
E„(V) has a resonance at coq and may have other 


co2-n2fi2 


(13b) 


where we have used the same notation as in equations 
(10) and (11). Here has an extra pole at U (and 
subharmonics). Despite this limitation, the forced 
oscillator is attractive because of its experimental 
simplicity. The choice of a particular mode of operation 
depends on the system under consideration, the 
requirements, and the instrumentation available. 


Experimental 

In keeping with former sections of this report, where 
we have discussed various possibilities for matching with 
specific situations, we shall outline here an experimental 
system flexible enough to be used in various modes of 
operation. 

The main requirements from the magnetic field B are 
intensity and homogeneity in the region of oscillations. A 
strong magnetic field creates high fundamental frequency 
for a given magnetic moment and inertia and thus makes 
possible a wider range of velocities. Also, the stronger the 
field, the higher the pickup amplitude and consequently 
the signal-to-noise ratio. For a cylindrical oscillator of 
1-g cm2 moment of inertia, a magnetic field of 104 q, and 
a magnetic moment of 1 G (a coil of l-cm2 cross-sectional 
area and 10-A turns), the fundamental resonance 
frequency wo/27r is about 17 Hz. This frequency, which 
corresponds to a surface peak tangential velocity of the 
order of 50 cm sec’-i, is achievable with quite moderate 
magnetic fields and moments. Enhancing the frequency is 
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not technically difficult, but one should be careful to 
keep the inertia term I6/R oc l/R o)2 smaller than the 
sliding friction force; otherwise, slip may occur. 

Even though the sliding friction force puts an upper 
limit on the fundamental frequency of oscillation, the 
analysis of higher harmonics is limited only by the rate of 
data acquisition and the time constant of the electronics. 
With modern equipment a data acquisition rate of 10-3 
relative resolution for both the oscillation period and 
voltage amplitude should not be difficult to achieve. 
Thus with careful shielding and grounding, very high 
order harmonics are resolvable. This feasible harmonics 
resolution provides valuable information on intrinsic 
frictional relaxation times much shorter than the 
oscillation period. Such may be, for example, the 
relaxation time needed for the friction to switch polarity 
around the zero velocity point. If this time is in the range 
10 - " sec, it could be detected by relatively standard 
equipment. To go to the other extreme of long relaxation 
times, the fundamental frequency should be lowered (to 
the 0.1 -Hz range). 

With approximately the magnetic field and the 
frequency mentioned above, the pickup coil can be 
designed to yield a signal amplitude of a few millivolts 
without exceeding the time-constant limit. This would 
allow a signal resolution better then 10-3. Altogether, 
instrumentation is not expected to be the cause of 
substantial experimental errors. All of the other errors 
foreseen can be classified into two groups: (1) errors 
stemming from gravitational-geometrical origins such as 
a possible lack of balance of either the oscillator or the 
planar surface and (2) errors resulting from uncertainty 
in determining the constants of the problem (the a/’s). 
The first group is not hard to reduce, by fine machining 
and mechanical balancing, to a minute value. Also, we 
can show that this type of error, when small, results in a 
small constant additive to gq. Thus it translates to be an 
uncertainty of the second group but a very small one. The 
additive constant would not cause any additional error to 
that which exists anyway in measuring gq (or coq). 

The other errors result from uncertainties in the values 
of the coefficients uq, ui, and G 2 (or coq» ^i)- For 

all nonconducting materials aj=0 (or can be directly 
measured if caused by feedback). For conductors the best 
policy would be to start with minimal induced currents. 
This objective could be achieved if the bulk of the 
oscillator is made of an insulator, where the tested 
sample, as a thin foil of thickness d, is attached to the 
angular section that makes contact during oscillations. 
With this geometry the induced current is minimized to 
the order of dB{Rd^)2aB0/c and the resulting force is 
{dAa) {B/c)W66^R, where A is the foiFs area, IRO^^j is 
the foil’s width, and a is the conductivity. The last force 
is somewhat different in character from the “Newtonian 


viscous” term discussed so far, which approximates the 
homogeneous bulk conductor, or the full cylindrical foil. 
Inserting values as in the former example and a =104 
(fl cm) - 1 and = 10 - 2 cm in the last expression, we get a 
force of 10-5 to 10-6 g. This value is negligible when 
compared with any typical rolling friction force, and it 
would not become significant even for the lowest 
resistivity metal. Only when the whole oscillator is made 
of a good conductor might the induced current produce a 
significant error, and then it should be treated along the 
lines described in the section Analysis. 

Since in most practical cases we can reduce Uj to a 
negligible value, there are two independent parameters to 
measure. One feasible possibility is to measure wo and gq. 
Note that Gq = -N^-A^If/NpAp, where and A^. are the 
number of turns and cross-sectional area of the restoring 
coil and is the current through the restoring coil. Thus 
Go is proportional to the ratio of the pickup voltages of 
these two coils when the cylinder is revolving at a 
constant angular velocity. This quantity can be measured 
to a very good accuracy. The parameter wq can be 
approximated by the actual measured frequency. If this is 
not accurate enough, a correction could be made by 
taking into account the loss of the first interation and 
making a second iteration calculation. To conclude, we 
estimate that absolute accuracies of the order of ±1 
percent can be obtained. For a comparison between two 
samples, the relative accuracy could be higher by one 
order of magnitude. 

A block diagram of the experimental system is shown 
in figure 2. The system makes possible the choice between 
feedback-maintained steady-state oscillations, pulse- 
excited decaying oscillations, or continuous-wave forced 
oscillations. The feedback control unit might include a 
filter if one chooses to feed back only the first harmonic. 

The oscillating body cannot be a sphere because any 
rotational component around the axis of the pickup coil 
is forbidden. There should be a cylindrical symmetry 
around the axis mutually normal to the coil’s axes. Apart 
from this, a great deal of freedom in the oscillator shape 
is allowed. It may be of any shape that has a sectional 
cylindrical symmetry. A measurement on a rolling bail 
could be done if only one degree of freedom is eliminated 
(e.g., by a rigid attachment of two identical balls). The 
finishing of the surfaces should be accurate enough to 
avoid microscopic slip. If the friction versus load has to 
be measured, we recommend varying the load by 
attaching weights on the oscillator. Loading by means of 
any mechanical system is undesirable since it may 
introduce tangential forces additional to those described 
in equation (3). Finally, since no assumption has been 
made on /, the experiment can be performed in any 
medium or in vacuum and on either dry or lubricated 
surfaces. 
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Oscilloscope 



Figure 2, - Block diagram of experimental system. 


Concluding Remarks 

The method suggested herein for measuring rolling 
friction is characterized by mechanical simplicity and is 
compatible with friction measurements and detecting 
relaxation times in the range lOi to 10-4 sec. It is, 
however, a relatively low- velocity measurement because 
of possible slip at higher velocities. We believe it is 
capable of detecting friction polarity switching times and 
friction variations that may occur in the low speed range 
(0.1 cm/sec, ref. 2). Further analysis is needed for 
interpreting the output spectral response of rolling 
friction torque in terms of constant -velocity friction. At 
low vibrational frequencies, when the acceleration is 
small enough not to smear out friction velocity- 
dependence variations, translation to velocity depen- 
dency is easy. In this case we can write 

/W=— V D„s\nno>t{e) (14) 

’T n=0,S... 
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7T 



COS 


where / ( 0) = r ( F) is the inverse of K( /) and t varies over 
one oscillation period. If high acceleration results in f{6) 
different from the steady-state function, we can still learn 
a great deal about the relaxation time required to 
establish the steady-state function. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 15, 1983 
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